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We investigate the fusion barriers for reactions involving Ca isotopes 40 Ca + 40 Ca, 40 Ca + 48 Ca, 
and 48 Ca + 48 Ca using the microscopic time-dependent Hartree-Fock theory coupled with a density 
constraint. In this formalism the fusion barriers are directly obtained from TDHF dynamics. We 
also study the excitation of the pre-equilibrium GDR for the 40 Ca + 48 Ca system and the associated 
7-ray emission spectrum. Fusion cross-sections are calculated using the incoming-wave boundary 
condition approach. We examine the dependence of fusion barriers on collision energy as well as on 
the different parametrizations of the Skyrme interaction. 
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I. INTRODUCTION 

The microscopic study of nuclear many-body problem 
and the understanding of the nuclear interactions that 
reproduce the observed structure and reaction proper- 
ties are the underlying challenges of low energy nuclear 
physics. In this context, detailed investigations of the 
fusion process will lead to a better understanding of the 
interplay among the strong, Coulomb, and weak interac- 
tions as well as the enhanced correlations present in these 
many-body systems. 

Recently, particular experimental attention has been 
given to fusion reactions involving Ca isotopes [H-Q. 
These new experiments supplemented the older fusion 
data [f| and extended it to lower sub-barrier energies. 
Comparison of the sub-barrier cross-sections with those 
calculated using standard coupled-channel calculations 
suggested a hindrance of the fusion cross-sections at deep 
sub-barrier energies [IH3 ■ One of the underlying reasons 
for the failure of standard coupled-channel approach is 
the use of frozen densities in the calculation of double- 
folding potentials, resulting in potentials that behave in 
a completely unphysical manner for deep sub-barrier en- 
ergies. While the outer part of the barrier is largely de- 
termined by the early entrance channel properties of the 
collision, the inner part of the potential barrier is strongly 
sensitive to dynamical effects such as particle transfer 
and neck formation. This has been remedied in part by 
extensions of the coupled-channel approach to include a 
repulsive core [6j or the incorporation of neck de gre es of 
freedom @, ||| . More recent calculations [3, d, [l(| us- 
ing the coupled-channel approach with a repulsive core 
have provided much improved fits to the data. A detailed 
microscopic study of the fusion process for Ca based re- 
actions 40 Ca + 40 Ca, 40 Ca + 48 Ca, and 48 Ca + 48 Ca could 
provide further insight into the reaction dynamics as well 
as a good testing ground for the theory since these iso- 
topes are commonly used in fitting the parameters of the 
effective nuclear interactions, such as the Skyrme force. 

During the past several years, we have developed a mi- 
croscopic approach for calculating heavy-ion interaction 



potentials that incorporates all of the dynamical entrance 
channel effects included in the time-dependent Hartree- 
Fock (TDHF) description of the collision process 
The method is based on the TDHF evolution of the nu- 
clear system coupled with density-constrained Hartree- 
Fock calculations (DC-TDHF) to obtain the ion-ion in- 
teraction potential. The formalism was applied to s tudy 
fusion cross-sections for the systems 132 Sn+ 64 Ni (Ra j. 

64 Ni+ 64 N j FJ| 16 O + 208p b ^ 132,124 Sn+ 96 Zr ^ ag 

well as to the study of the entrance channel dynamics 
of hot and cold fusion reactions leading to superheavy 
element Z — 112 fl6| . and dynamical excitation ener- 
gies In all cases, we have found good agreement 
between the measured fusion cross sections and the DC- 
TDHF results. This is rather remarkable given the fact 
that the only input in DC-TDHF is the Skyrme effective 
N-N interaction, and there are no adjustable parameters. 

In Section[TT]we outline the main features of our micro- 
scopic approach, the DC-TDHF method. In Section |TT] 
we also discuss the calculation of ion-ion separation dis- 
tance, coordinate-dependent mass, calculation of fusion 
cross-sections, and giant dipole resonance (GDR) formal- 
ism. In Sec. IIIII we present interesting aspects of the 
reaction dynamics and compare our results with exper- 
iment and other calculations. In Sec. [V] we summarize 
our conclusions. 



II. FORMALISM 

A. DC-TDHF Method 

In the DC-TDHF approach [11] the TDHF time- 
evolution takes place with no restrictions. At certain 
times during the evolution the instantaneous density 
is used to perform a static Hartree-Fock minimization 
while holding the neutron and proton densities con- 
strained to be the corresponding instantaneous TDHF 
densities [l||, Il9| . In essence, this provides us with 
the TDHF dynamical path in relation to the multi- 
dimensional static energy surface of the combined nu- 
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clear system. The advantages of this method in com- 
parison to other mean-field based microscopic methods 
such as the constrained Hartree-Fock (CHF) method are 
obvious. First, there is no need to introduce artificial 
constraining operators which assume that the collective 
motion is confined to the constrained phase space: sec- 
ond, the static adiabatic approximation is replaced by 
the dynamical analogue where the most energetically fa- 
vorable state is obtained by including sudden rearrange- 
ments and the dynamical system does not have to move 
along the valley of the potential energy surface. In short 
we have a self-organizing system which selects its evolu- 
tionary path by itself following the microscopic dynam- 
ics. All of the dynamical features included in TDHF are 
naturally included in the DC-TDHF calculations. These 
effects include neck formation, mass exchange, internal 
excitations, deformation effects to all order, as well as 
the effect of nuclear alignment for deformed systems. In 
the DC-TDHF method the ion-ion interaction potential 
is given by 

V(R) = Et> c (R)-E Ai -E A2 , (1) 

where -Edc is the density-constrained energy at the in- 
stantaneous separation R(t), while E Al and Ea 2 are the 
binding energies of the two nuclei obtained with the same 
effective interaction. In writing Eq. (TTJ) we have intro- 
duced the concept of an adiabatic reference state for a 
given TDHF state. The difference between these two en- 
ergies represents the internal energy. The adiabatic refer- 
ence state is the one obtained via the density constraint 
calculation, which is the Slater determinant with lowest 
energy for the given density with vanishing current and 
approximates the collective potential energy 18]. We 
would like to emphasize again that this procedure does 
not affect the TDHF time-evolution and contains no free 
parameters or normalization. 

In addition to the ion-ion potential it is also possible 
to obtain coordinate dependent mass parameters. One 
can compute the "effective mass" M(R) using the con- 
servation of energy 

M («, = *isy>i , p, 

where the collective velocity R is directly obtained from 
the TDHF evolution and the potential V(R) from the 
density constraint calculations. In calculating fusion 
cross-sections this coordinate-dependent mass is used to 
obtain a transformed ion-ion potential as described be- 
low. 



B. Excitation Energy 

The calculation of the excitation energy is achieved by 
dividing the TDHF motion into a collective and intrin- 
sic part. The major assumption in achieving this goal 



is to assume that the collective part is primarily deter- 
mined by the density p(r, t) and the current j(r, t). Con- 
sequently, the excitation energy can be formally written 
as [13 

E*(t) = E TDHF - E cM (p(t)J(t)) , (3) 

where -Etdhf is the total energy of the dynamical system, 
which is a conserved quantity, and £coii represents the 
collective energy of the system. In the next step we break 
up the collective energy into two parts 

E coU (t) = E kin (p(t), j(t)) + E DC {pit)) , (4) 

where -Ekin represents the kinetic part and is given by 

E kin (Pit), j(t)) = f / d 3 rf(t)/p(t), (5) 

which is asymptotically equivalent to the kinetic energy 
of the relative motion, ^pR 2 , where p is the reduced 
mass and R(t) is the ion-ion separation distance. The 
dynamics of the ion-ion separation R(t) is provided by an 
unrestricted TDHF run thus allowing us to deduce the 
excitation energy as a function of the distance parameter, 
E*(R). 

C. Calculation of R 

In practice, TDHF runs are initialized with energies 
above the Coulomb barrier at some large but finite sepa- 
ration. The two ions are boosted with velocities obtained 
by assuming that the two nuclei arrive at this initial sep- 
aration on a Coulomb trajectory. Initially the nuclei are 
placed such that the point x — in the x — z plane 
is the center of mass. During the TDHF dynamics the 
ion-ion separation distance is obtained by constructing a 
dividing plane between the two centers and calculating 
the center of the densities on the left and right halves of 
this dividing plane. The coordinate R is the difference 
between the two centers. The dividing plane is deter- 
mined by finding the point at which the tails of the two 
densities intersect each other along the x-axis. Since the 
actual mesh used in the TDHF calculations is relatively 
coarse we use a cubic-spline interpolation to interpolate 
the profile in the x-direction and search for a more pre- 
cise intersection value. This procedure has been recently 
described in Ref. [13] in great detail. 

The standard procedure for calculating R as described 
above starts to fail after a substantial overlap is reached 
(for R values smaller than the ones studied in this 
manuscript). We have observed that if one defines the 
ion- ion separation as R = Roy\Qw\, where Q20 is the 
mass quadrupole moment for the entire system, calcu- 
lated by using the collision axis as the symmetry axis, 
and Rq is a scale factor determined to give the correct 
initial separation distance at the start of the calculations. 
Calculating R this way yields almost identical results to 
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the previous procedure until that procedure begins to fail 
and continues smoothly after that point. Of course the 
minimum value of R calculated this way is never zero but 
is determined by the quadrupole moment of the compos- 
ite system. 



D. Fusion Cross-Section 

We now outline the calculation of the total fusion cross 
section using an arbitrary coordinate-dependent mass 
M(R). Starting from the classical Lagrange function 
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L{R,R) = ^M(R)R 2 



V(R), 

we obtain the corresponding Hamilton function 

P 2 



H(R,P) 



2M(R) 



V(R) 



(6) 



(7) 



where the canonical momentum is given by P = M(R)R. 
Following the standard quantization procedure for the 
kinetic energy in curvilinear coordinates [2lj 



T : 



2 



fJ 



d 



d 

dq v 



(8) 



where g^ v {q) denotes the metric tensor and g^ v {q) the 
reciprocal tensor, one obtains the quantized Hamiltonian 



H(R,P) = - M{R)-*PM{R)-*P 



with the momentum operator P 
fusion cross cross-section 



V(R) . (9) 
ihd/dR. The total 



As a result of this point transformation, both the clas- 
sical Hamilton function, Eq. ([7]), and the corresponding 
quantum mechanical Hamiltonian, Eq. ([9]), now assume 
the form 



P 2 



U{R) , 



(14) 



and the scaled heavy-ion potential is given by the expres- 
sion 

U{R) = V{R) = V{f- 1 {R)). (15) 

The fusion barrier penetrabilities T^(i? c . m .) are ob- 
tained by numerical integration of the two-body 
Schrodinger equation 



-ft 2 d 2 L{L + l)h 2 



2/z dR? 



2/ii? 2 



U{R) 



i> = 0, (16) 



using the incoming wave boundary condition (IWBC) 
method (23|. IWBC assumes that once the minimum 
of the potential is reached fusion will occur. In prac- 
tice, the Schrodinger equation is integrated from the po- 
tential minimum, R m i n , where only an incoming wave 
is assumed, to a large asymptotic distance, where it 
is matched to incoming and outgoing Coulomb wave- 
functions. The barrier penetration factor, T^{E c m ) is 
the ratio of the incoming flux at i? m i n to the incoming 
Coulomb flux at large distance. Here, we implement the 
IWBC method exactly as it is formulated for the coupled- 
channel code CCFULL described in Ref. [3] ■ This gives 
us a consistent way for calculating cross-sections at above 
and below the barrier energies. 



E. GDR Excitation 
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L=0 



l)T L 



(10) 



can be obtained by calculating the potential barrier pen- 
etrabilities Tl from the Schrodinger equation for the rela- 
tive motion coordinate R using the Hamiltonian ([9]) with 
an additional centrifugal potential 



H(R,P) 



H 2 L{L + 1) 
2M(R)R 2 



-E c 



il> L {R) = . (11) 



Alternatively, instead of solving the Schrodinger equa- 
tion with coordinate dependent mass parameter M(R) 
for the heavy-ion potential V(R), we can instead use the 
constant reduced mass fi and transfer the coordinate- 
dependence of the mass to a scaled potential U (R) using 
the well known coordinate scale transformation [22J 

dH=(^)*dfl. (12) 

Integration of Eq. (fl2|) yields 

R = f(R) R = f- 1 (R). (13) 



Let us now consider a central collision in x-direction 
and introduce the quantity 

1 Z 1 N 



D(t) 



NZ 



P =i 



71 = 1 



(17) 

which represents the expectation value of the as- 
component of the dipolc operator d x /e taken with the 
time-dependent TDHF Slater determinant |4>(f) >. Fol- 
lowin g th e bremsstrahlung approach developed by Baran 
et al. [25l [26j j we define the dipole acceleration 



D"{t) = 



d 2 D{t) 
dt 2 



(18) 



(19) 



and introduce its Fourier transform 

D"{uj) = J™* D"(t)e iut dt. 

Alternatively, for nearly harmonic vibrations one can use 
the expression \D"(uj)\ 2 = lo a \D(lu)\ 2 , with 

D(lj) = / tmaX D(t)e^ sin 4 ( n * ~ Wn ) dt . (20) 
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The time filtering sin is used to smooth out peaks com- 
ing from finite integration time. The "power spectrum" 
of the electric dipole radiation is given by (2f| 



60 



dP 



2a 



(21) 



where a = e 2 /(he) k, 1/137 denotes the fine structure 
constant. Recently, pre-equilibrium GDR excitation has 
also been studied in the context of TDHF 123 • 



III. 



RESULTS 



Calculations were done in 3-D geometry and using 
the full Skyrme interaction including all of the time- 
odd terms in the mean-field Hamiltonian (28[. The pri- 
mary Skyrme parametrization used was SLy4 |29j but we 
have also tested the new UNEDFO [13] and UNEDF1 [3l| 
parametrizations. For the reactions studied here, the lat- 
tice spans 48 fm along the collision axis and 15 fm in the 
other two directions. Derivative operators on the lattice 
are represented by the Basis-Spline collocation method. 
One of the major advantages of this method is that we 
may use a relatively large grid spacing of 1.0 fm and nev- 
ertheless achieve high numerical accuracy. The initial 
separation of the two nuclei is 18 fm for central colli- 
sions. The time-propagation is carried out using a Tay- 
lor series expansion (up to orders 10 — 12) of the unitary 
mean-field propagator, with a time-step At = 0.4 fm/c. 
We have performed density constraint calculations every 
10 — 20 time steps. The accuracy of the density con- 
straint calculations is commensurate with the accuracy 
of the static calculations. 



A. 40 Ca+ 40 Ca System 

In this subsection we will present our results for the 
40 Ca+ 40 Ca system. Most of our general discussions will 
be provided here as they would be the same for other 
systems. Specific points about individual systems and 
comparison of results for the three systems will be taken 
up in the subsequent sections. 

In Fig. Q] we show the microscopic ion-ion potential bar- 
riers obtained using Eq. (JTJ) calculated at three different 
collision energies, Stdhf = 55 MeV (black solid curve), 
-Etdhf = 60 MeV (red solid curve), and -Etdhf = 
65 MeV (blue solid curve). We observe a relatively small 
dependence on the collision energy. At lower energies 
the system has more time available for rearrangements to 
take place through the formation of a neck, whereas this 
is less and less the case at higher energies and the poten- 
tial barrier approaches the frozen-density limit [20] . This 
produces the observed trend in Fig. [TJ where the lowest 
barrier peak corresponds to the lowest energy and the 
barrier height increases with increasing collision energy. 
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FIG. 1. (Color online) Potential barriers, V(R), for the 
40 Ca-|- 40 Ca system obtained from density constrained TDHF 
calculations using Eq. (ff]) at three different energies, -Btdhf = 
55 MeV (black solid curve), £ T dhf = 60 MeV (red solid 
curve), and -Etdhf = 65 MeV (blue solid curve). The three 
dashed curves correspond to the transformed potential in 
Eq. (|15[) using the coordinate dependent masses. Also shown 
is the point Coulomb potential. 



The barrier heights are 53.02, 53.43, and 53.57 MeV, re- 
spectively. Similar energy dependence was also observed 
in the DD-TDHF calculations of Ref. [2(| and becomes 
more prevalent for heavier systems. Similarly, the po- 
sition of the barrier peak moves toward smaller R val- 
ues, albeit very slowly in this case, with increasing en- 
ergy. Corresponding R values for the barrier maximum 
are 10.41, 10.32, and 10.23 fm. What is also shown on 
Fig. [1] is the Coulomb potential assuming the two nuclei 
to be point particles with Z = 20. During the approach 
phase the microscopically calculated DC-TDHF poten- 
tial traces the point Coulomb potential, differing by less 
than 150 keV, which provides a test for the numerical 
accuracy. 

The energy dependence of potential barriers can also 
be understood if we examine the coordinate dependent 
mass of Eq. |(3J) shown in Fig. [5] for three different en- 
ergies, -Etdhf = 55 MeV (black solid curve), -Etdhf = 
60 MeV (red solid curve), and £ T dhf = 65 MeV (blue 
solid curve). The independence of this mass at lower 
energies is very similar to the one found in CHF calcula- 
tions (22|. On the other hand, at higher energies the co- 
ordinate dependent mass essentially becomes flat, which 
is again a sign that most dynamical effects are contained 
at lower energies. The peak at small R values is due to 
the fact that the center-of-mass energy is above the bar- 
rier and the denominator of Eq. ([2} becomes small due to 
the slowdown of the ions. We have used the coordinate 
dependent masses shown in Fig. [5] to obtain the scaled 
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FIG. 2. (Color online) Coordinate dependent mass M(R) 
scaled by the constant reduced mass fi, obtained from Eq. 
at three different TDHF energies. 



are shown in Fig. [3] We observe that all of the scaled 
barriers give a very good description of the experimental 
fusion cross-sections. The high energy part of the fu- 
sion cross-sections are primarily determined by the bar- 
rier properties in the vicinity of the barrier peak. On 
the other hand sub-barrier cross-sections are influenced 
by what happens in the inner part of the barrier and 
here the dynamics and consequently the coordinate de- 
pendent mass becomes very important. As we observe in 
Fig.[3]we get a very good agreement with experiment for 
the barriers obtained at the lowest two collision energies, 
whereas the E c , m , — 65 MeV curve slightly underesti- 
mates the cross-section at lower energies. Although not 
shown in Fig. [3] the cross-sections obtained using the the 
unsealed potentials V(R) and a constant reduced mass 
/i also agree well with the data at higher energies but 
either significantly over-estimate or under-estimate the 
cross-section at lower energies due to the absence of the 
coordinate dependent mass. In summary, the calculated 
fusion cross-sections for the 40 Ca+ 40 Ca system reproduce 
the experimental cross-sections reasonably well, which is 
a testament that TDHF with Skyrme force provides a 
good description for this collision. 



potentials U(R) of Eq. (|15p . These potentials are shown 
as the dashed curves in Fig. [TJ As we see the coordinate 
dependent mass only changes the inner parts of the bar- 
riers for all energies. Furthermore, the effect is largest for 
the lowest energy collision and diminishes as we increase 
the collision energy. 
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FIG. 3. (Color online) Total fusion cross section as a function 
of -Ec.m.- Three separate theoretical cross section calculations 
are shown, based on the energy-dependent DC-TDHF heavy- 
ion potentials V(R) at energies Etdhf = 55, 60, 65 MeV. 
The experimental data (filled circles) are taken from Ref. Q . 
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FIG. 4. (Color online) Potential barriers, V(R), for the 
48 Ca+ 48 Ca system obtained from density constrained TDHF 
calculations using Eq. (JXJ) at three different energies, Etdhf — 
55 MeV (black solid curve), -Etdhf = 65 MeV (red solid 
curve), and -Etdhf = 70 MeV (blue solid curve). The dashed 
curves correspond to the transformed potential in Eq. (|15[) 
using the coordinate dependent masses. Also shown is the 
point Coulomb potential. 



B. 48 Ca+ 48 Ca System 



We have obtained the fusion cross sections by numer- 
ical integration of Eq. (flU)) . The resulting cross-sections 



In Fig. d] we show the microscopic ion-ion potential 
barriers for the stable neutron- rich 48 Ca+ 48 Ca system 
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calculated at three different collision energies, -Etdhf = 
55 MeV (black solid curve), E TDHF = 65 MeV (red solid 
curve), and -Etdhf = 70 MeV (blue solid curve). The 
barrier heights, for increasing collision energy, are 50.98, 
51.37, and 51.52 MeV, respectively, with corresponding R 
values of 10.80, 10.63, and 10.52 fm. The comparison of 
these barriers with those of the 40 Ca+ 40 Ca system shows 
that the barrier heights are reduced by about 2 MeV and 
the location of the barrier maximum is at a slightly larger 
R value. This is due to the fact that two two 48 Ca nuclei 
are larger than the corresponding 40 Ca nuclei and thus 
their outer skins come into contact at a larger R value. 
After this point the nuclear interaction sets in causing 
the trajectory to deviate from the point Coulomb one 
and producing a peak at a lower energy value. 
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FIG. 5. (Color online) Fusion cross sections for the 
48 Ca+ 48 Ca system as a function of E c . m .. Three separate 
theoretical cross section calculations are shown, based on the 
energy-dependent DC-TDHF heavy-ion potentials V(R) at 
energies -Etdhf = 55, 65, and 70 MeV. The experimental 
data (filled circles) are taken from Ref. 0]. 

Figure [5] shows the fusion cross-sections for the 
48 Ca+ 48 Ca system calculated using the barriers of Fig.|4] 
While the overall quality of the agreement with the ex- 
perimental fusion cross-sections is good, specially for the 
last two collision energies, the curve corresponding to 
the lowest collision energy of 55 MeV overestimates the 
experiment for lower -E cm . values. This illustrates the 
sensitivity of the results to the height of the potential 
barrier, the difference in this case between the two bar- 
riers being around 0.4 MeV. 

While the quality of the DC-TDHF results for the 
48 Ca+ 48 Ca system is very good for a parameter free 
microscopic approach, we have decided to investigate 
this further. One of the problems is that the Skyrme 
fits do badly in reproducing the single-particle proper- 
ties [32J, particularly the neutron single-particle states. 
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FIG. 6. (Color online) Total fusion cross section as a func- 
tion of -Bc.m. for two different parametrizations of the Skyrme 
force, SLy4 and UNEDF1, at the collision energy of -Etdhf = 
55. The experimental data (filled circles) are taken from 
Ref. ||. 

For example, with SLy4 parametrization the neutron lev- 
els ld 5 / 2 , 2si/2, and ld 3 / 2 in 48 Ca nucleus are about 
7 — 3 MeV lower in energy than the corresponding ex- 
perimental values [3p| . In one of the parametrizations of 
the nuclear density functional, UNEDF0, these single- 
particle energies are raised to values closer to experi- 
mental ones. However, in this case the magic gap at 
N = 28 vanishes [3(| ■ Recently, a new parametrization, 
UNEDF1, was introduced with a better incorporation of 
large deformations among other improvements |3l| . An- 
other possible advantage of the UNEDF1 parametriza- 
tion for TDHF calculations is that no center-of-mass cor- 
rection term was used in the functional, which is also 
the case in TDHF. We have tried this parametrization 
in our DC-TDHF calculation for the 48 Ca+ 48 Ca system 
at E c .m. = 55 MeV. In practice, we had to reduce our 
time-step from At = 0.4 fm/c to At = 0.1 fm/c, and our 
density-constraint convergence parameters by a factor of 
ten or more. This is probably due to the larger power 
of the density in the £3 term of the Skyrme interaction. 
The peak of the fusion barrier is slightly broader and 
higher in energy by 187 keV for UNEDF1. The use of 
UNEDF1 results in an improvement of the fusion cross- 
sections as shown in Fig. [5] We have also tried using the 
parametrization SLy5 but this resulted in no appreciable 
difference from the SLy4 case. 



C. 40 Ca+ 48 Ca System 

In this section we will examine the fusion of the asym- 
metric 40 Ca+ 48 Ca system. For highlighting the differ- 
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FIG. 7. (Color online) Experimental fusion cross sections as 
a function of E c . m . for the three systems studied. The exper- 
imental data are taken from Ref. [J]. 



ences among all three systems we have plotted only the 
experimental fusion cross-sections in Fig. [71 Relative to 
the 40 Ca+ 40 Ca and 48 Ca+ 48 Ca systems the 40 Ca+ 48 Ca 
fusion cross-sections show a different systematic behav- 
ior. At higher bombarding energies the 40 Ca+ 48 Ca cross- 
sections fall below the 40 Ca-|- 4cl Ca data points, whereas 
at sub-barrier energies they rise above the 48 Ca+ 48 Ca 
data. Historically, the theoretical description of the fu- 
sion cross-sections for the 40 Ca+ 48 Ca system have been 
complicated by the couplings to various transfer channels 
with positive Q values [33l . l34| . A strong enhancement of 
fusion cross-sections at lower energies seen in Fig. [JJ was 
attributed to this effect 0. 

In Fig. [5] we show the microscopic ion-ion potential 
barriers for the 40 Ca-|- 48 Ca system calculated at three 
different collision energies, £tdhf = 55 MeV (black 
solid curve), -Etdhf = 60 MeV (red solid curve), and 
-Etdhf = 70 MeV (blue solid curve). The barrier 
heights, for increasing collision energy, are 51.11, 51.45, 
and 51.81 MeV, respectively, with corresponding R val- 
ues of 10.75, 10.65, and 10.51 fm. We have used the 
coordinate dependent masses to obtain the scaled poten- 
tials U (R) of Eq. (TTBl . These potentials are also shown 
as the dashed curves in Fig. [5J As before the coordi- 
nate dependent mass only changes the inner parts of the 
barriers for all energies and the effect diminishes as we 
increase the collision energy. 

Figure |H] shows the fusion cross-sections for the 
40 Ca+ 48 Ca system calculated using the barriers of Fig. [5] 
The observed trend for sub-barrier energies is typical for 
DC-TDHF calculations when the underlying microscopic 
interaction gives a good representation of the participat- 
ing nuclei. Namely, the potential barrier corresponding 
to the lowest collision energy gives the best fit to the 
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FIG. 8. (Color online) Potential barriers for the 40 Ca+ 48 Ca 
system obtained from density constrained TDHF calculations 
at three different energies, Etdhf = 55 MeV (black solid 
curve), -Etdhf = 60 MeV (red solid curve), and -Etdhf = 
70 MeV (blue solid curve). The three dashed curves corre- 
spond to the transformed potential in Eq. (|15l) using the co- 
ordinate dependent masses. Also shown is the point Coulomb 
potential. 
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FIG. 9. (Color online) Fusion cross sections for the 
40 Ca+ 48 Ca system as a function of E c . m .- Three separate 
theoretical cross section calculations are shown, based on the 
energy-dependent DC-TDHF heavy-ion potentials V(R) at 
energies Etdhf = 55, 60, and 70 MeV. The experimental 
data (filled circles) are taken from Ref. 0]. 



sub-barrier cross-sections since this is the one that al- 
lows for more rearrangements to take place and grows the 
inner part of the barrier. Considering the fact that his- 
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torically the low-energy sub-barrier cross-sections of the 
40 Ca+ 48 Ca system have been the ones not reproduced 
well by the standard models, the DC-TDHF results are 
quite satisfactory, indicating that the dynamical evolu- 
tion of the nuclear density in TDHF gives a good overall 
description of the collision process. The shift of the cross- 
section curve with increasing collision energy is typical. 
In principle one could perform a DC-TDHF calculation at 
each energy above the barrier and use that cross-section 
for that energy. However, this would make the compu- 
tations extremely time consuming and may not provide 
much more insight. 

The trend at higher energies is atypical. The cal- 
culated cross-sections are larger than the experimental 
ones by about a factor of two. Such lowering of fusion 
cross-sections with increasing collision energy is com- 
monly seen in lighter systems where various inelastic 
channels, clustering, and molecular formations are be- 
lieved to be the contributing factors |35|. In the recent 
coupled-channel approach this is solved by the addition 
of a small ima gina ry potential near the minimum of the 
repulsive core [36| . Such an imaginary part was found 
not to be necessary for the 40 Ca+ 40 Ca and 48 Ca+ 48 Ca 
systems d, ■ At this time we do not have access to 
coupled-channel results for the 40 Ca+ 48 Ca system. We 
have repeated our DC-TDHF calculations using the UN- 
EDF1 interaction which resulted in a small improvement, 
reducing the difference with the experimental values to 
a factor of about 1.5. This issue will be discussed fur- 
ther in the next Section where we examine the excitation 
properties obtained from TDHF calculations. 
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FIG. 10. (Color online) Neutron and proton transfer as a 
function of ion- ion separation distance R for the 40 Ca+ 4S Ca 
system. Solid lines denote the particles originally belonging 
to the 48 Ca nucleus and dashed lines to the 40 Ca nucleus. The 
dotted vertical line shows the location of the potential barrier 
peak, Rb- 



of the system from quasi-fission. The dynamical excita- 
tion energies calculated from DC-TDHF are relative to 
an unequilibrated composite system rather than a true 
compound nucleus. 



IV. EXCITATIONS 

Excitations are believed to have a significant impact 
on the outcome of the fusion reactions. The excitations 
can range from the entrance channel quantal excitations 
of the projectile and target, as in the coupled-channel 
approach, to collective excitations of pre-equilibrium sys- 
tem, to compound nucleus excitations. These can be fur- 
ther influenced by particle transfer, pre-equilibrium emis- 
sions, and evaporation, among others. Theoretically such 
effects are commonly introduced by hand into various re- 
action models. However, the influence of excitations on 
nuclear reaction dynamics remains to be a difficult and 
an open problem as it combines both nuclear structure 
and dynamics under nonequilibrium conditions. 

In Sec. Ill Bl we have outlined the calculation of the 
dynamical excitation energy within the DC-TDHF for- 
malism. This approach was used to calculate excitation 
energies for various systems, including the excitation en- 
ergy of the heavy systems leading to superheavy forma- 
tions [ll|. Here, we have shown that what is also impor- 
tant is the excitation of the system at the point of cap- 
ture, which is a deciding factor for forming a composite 
system or fusion-fission. This is different than excitation 
of the compound nucleus, which determines the survival 
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FIG. 11. (Color online) Excitation energy, E*(R), as a func- 
tion of the ion-ion separation distance, R, for the three sys- 
tems studied here. 

The systems studied in this manuscript present inter- 
esting target-projectile combinations for studying excita- 
tions. For example, the non-symmetric 40 Ca+ 48 Ca sys- 
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tern will have a pre-equilibrium GDR excitation in com- 
parison to the other two symmetric systems, as well as 
particle transfer. Naturally, some of these effects are in- 
cluded in the fusion barriers obtained via the DC-TDHF 
method, as they influence the change in the TDHF den- 
sity used in the density constraint calculation. The quan- 
tity E*(R(t)) calculated from via TDHF and DC-TDHF 
represents the average of dynamical excitations present 
in the mean-field theory, with the exception of excita- 
tions not determined by the nuclear density and current 
(such as the spin currents). Another point to consider 
is the pairing interaction. It has been a general consen- 
sus that during a heavy-ion collision the effects of pairing 
are washed away due to the high excitations. In this par- 
ticular case paring effects are minimal for the initial Ca 
nuclei studied. However, even if the pairing can be ig- 
nored during the collision process, in many cases it plays 
an important role for obtaining good initial HF states 
as well as obtaining realistic density constraint solutions 
when a single composite is formed. Fortunately, the lat- 
ter case corresponds to the region of barrier minimum 
and not the region around the barrier peak where most 
fusion cross-sections are measured. This may in some 
way explain the success of the DC-TDHF barriers for 
fusion. 




1000 2000 3000 4000 5000 6000 7000 
t (fm/c) 

FIG. 12. (Color online) Time evolution of isoscalar deforma- 
tion parameter P2 for the head-on collision of all three systems 
at a collision energy of -Etdhf = 55 MeV (top panel), and the 
time evolution of the isovector dipole amplitude D(t) for the 
40 Ca+ 48 Ca system (bottom panel). 

In Fig.Hniwe plot the average number of neutrons and 
protons transferred during the early stages of the TDHF 
collision of the 40 Ca+ 48 Ca system at E TDHF = 60 MeV. 
The solid lines denote the neutrons and protons asymp- 
totically belonging to the 48 Ca nucleus and dashed lines 
to the 40 Ca nucleus. A number of interesting things can 
be observed from the plot; first is the fact that most of 



the transfer seems to start after we pass the potential bar- 
rier peak. This indicates that particle transfer primarily 
modifies the inner part of the barrier and not so much 
the barrier height. The other observation is that on av- 
erage about three neutrons are transferred from 48 Ca to 
40 Ca but there is also a small amount of proton transfer 
in the opposite direction. 
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FIG. 13. (Color online) Power spectrum of the isovector 
dipole amplitude D(t) for the 40 Ca-|- 48 Ca system is shown 
for various time-intervals 

Figure QT] shows the excitation energy, E*(R), for the 
three systems studies here. The excitation energy was 
calculated for the same value of e = E c , m ./fi — 2.75 MeV 
for all systems, which corresponds to collision energies 
of 55, 60, and 66 MeV, respectively. All curves initially 
behave in a similar manner, at large distances the exci- 
tation is zero, as the nuclei approach the barrier peak 
the excitations start and monotonically rise for larger 
overlaps. The interesting observation is that the excita- 
tions for the intermediary 40 Ca+ 48 Ca system start at a 
slightly earlier time and rise above the other two systems. 
This may be largely due to the fact that an asymmetric 
system has some additional modes of excitation in com- 
parison to the other two symmetric systems. It may be 
plausible to consider the direct influence of the excitation 
energy, E*(R), on the fusion barriers by making an anal- 
ogy with the coupled-channel approach and construct a 
new potential V*{R) = V{R) + E*(R), which has all 
the excitations added to the ion- ion potential V(R) that 
should be calculated at higher energies to minimize the 
nuclear rearrangements (frozen-density limit). The re- 
sulting potentials somewhat resemble the repulsive-core 
coupled-channel potentials of Ref. @- This approach 
does lead to improvements in cases where most of the 
excitation energy is in the form of collective excitations 
rather than irreversible stochastic dissipation (true espe- 
cially for lighter systems). For the cases studied here we 
do see a small improvement for the fusion cross-sections 



10 



of the Ca+ 48 Ca system. The viability of this approach 
requires further examination and will be studied in the 
future. 

One of the excitation modes for the 40 Ca+ 48 Ca sys- 
tem, namely the particle transfer, was already discussed 
above. The others are various isovector modes such as the 
pre-equilibrium GDR excitation. In Fig. [12] we show the 
time development of the isoscalar deformation parame- 
ter fa = ^f(r 2 Y 2 o}/Ar rms (top panel) and the isovec- 
tor dipole amplitude D(t) of Eq. (fTTj) (bottom panel). 
The isovector amplitudes for the two symmetric systems 
would appear as a zero-line on this plot. We have left 
out the initial approach phase of the collision (about 
260 fm/c) from these plots. The time-evolution of the 
deformation fa is very similar for all three systems ex- 
cept at larger times the 40 Ca+ 48 Ca system moves toward 
smaller deformation in comparison to the other two sys- 
tems. This is again most likely due to increased exci- 
tations that drive the system toward a more compact 
shape. The evolution of the dipole operator D(t) shows 
little damping over a long time interval. The power spec- 
trum associated with this time time-evolution is shown 
in Fig. [T31 We have evaluated the spectrum for differ- 
ent time-intervals. The broadest curve corresponds to a 
very early stage of the collision, 260 — 1000 fm/ c, but has 
the lowest amplitude (multiplied by 10 in Fig. [13]). In 
the second interval, 260 — 3000 fm/c, we see the develop- 
ment of various peaks and an increase in amplitude. The 
spectrum for the entire time interval, 260 — 7200 fm/c, 
is dominated by a sharp peak around 11 MeV, which is 
primarily originating from the later stages of the collision 
which is evident if one compares it to the spectrum from 



the interval 3000 - 7200 fm/c. 

V. SUMMARY 

In this manuscript we have provided a microscopic 
study of Ca+Ca fusion using the DC-TDHF approach. 
These reactions have recently been of considerable inter- 
est for the fusion community with a flurry of phenomeno- 
logical analyses of the fusion data. Here, we have pro- 
vided a microscopic alternative to this analysis. We have 
shown that microscopically obtained ion-ion potentials 
do give a reasonably good description of the fusion cross- 
sections. 

The fully microscopic TDHF theory has shown it- 
self to be rich in nuclear phenomena and continues to 
stimulate our understanding of nuclear dynamics. The 
time-dependent mean-field studies seem to show that the 
dynamic evolution builds up correlations that are not 
present in the static theory. While modern Skyrme forces 
provide a much better description of static nuclear prop- 
erties in comparison to the earlier parametrizations there 
is a need to obtain even better parametrizations that in- 
corporate deformation and reaction data into the fit pro- 
cess. 

ACKNOWLEDGMENTS 

This work has been supported by the U.S. Depart- 
ment of Energy under grant No. DE-FG02-96ER40963 
with Vanderbilt University. RK thankfully acknowl- 
edges The Scientific and Technological Council of Turkey 
(TUBITAK, BIDEB-2219) for a partial support. 



A. M. Stefanini et al, Phys. Lett. B 679, 95 (2009). 
C. L. Jiang et al, Phys. Rev. C 82, 041601 (2010). 
G. Montagnoli et al, Nucl. Phys. A 834, 159c (2010). 

G. Montagnoli and A. M. Stefanini, EPJ Web Conf. 17, 
05001 (2011). 

H. A. Aljuwair et al, Phys. Rev. C 30, 1223 (1984). 

S. Misicu and H. Esbensen, Phys. Rev. Lett. 96, 112701 
(2006). 

Takatoshi Ichikawa, Kouichi Hagino, and Akira Iwamoto, 
Phys. Rev. C 75, 057603 (2007). 

Takatoshi Ichikawa, Kouichi Hagino, and Akira Iwamoto, 
Phys. Rev. C 75, 064612 (2007). 
H. Esbensen et al, Phys. Rev. C 82, 054621 (2010). 
§erban Mi§icu and Florin Carstoiu, Phys. Rev. C 84, 
05160RR) (2011). 

Oberacker, 



E. 



Rev. C 74, 



A. S. Umar and V. 
021601(R) (2006). 

A. S. Umar and V. E. Oberacker, 
014614 (2007). 

A. S. Umar and V. E. Oberacker, 
064605 (2008). 

A. S. Umar and V. E. Oberacker, Eur. Phys. J. A 39, 
243 (2009). 

V. E. Oberacker, A. S. Umar, J. A. Maruhn, and P.-G. 



Phys. 
Phys. 
Phys. 



Rev. C 76, 



Rev. C 77, 



Reinhard, Phys. Rev. C 82, 034603 (2010). 

[16] A. S. Umar, V. E. Oberacker, J. A. Maruhn, and P.-G. 
Reinhard, Phys. Rev. C 81, 064607 (2010). 

[17] A. S. Umar, V. E. Oberacker, J. A. Maruhn, and P.-G. 
Reinhard, Phys. Rev. C 80, 041601(R) (2009). 

[18] R. Y. Cusson, P. -G. Reinhard, M. R. Strayer, J. A. 
Maruhn, and W. Greiner, Z. Phys. A 320, 475 (1985). 

[19] A. S. Umar, M. R. Strayer, R. Y. Cusson, P. -G. Rein- 
hard, and D. A. Bromley, Phys. Rev. C 32, 172 (1985). 

[20] Kouhei Washiyama and Denis Lacroix, Phys. Rev. C 78, 
024610 (2008). 

[21] MR. Spiegel, Vector Analysis, Schaum's Outline Series, 

McGraw-Hill (1959). 
[22] K. Goeke, F. Griimmer, and P. -G. Reinhard, Ann. Phys. 

150, 504 (1983). 
[23] G. H. Rawitscher, Phys. Rev. 135, 605 (1964). 
[24] K. Hagino, N. Rowley, and A. T. Kruppa, Comp. Phys. 

Comm. 123, 143 (1999). 
[25] V. Baran et al., Nucl. Phys. A 600, 111 (1996). 
[26] V. Baran et al., Phys. Rev. C 79, 021603(R) (2009). 
[27] C. Simenel, Ph. Chomaz, and G. de France, Phys. Rev. 

C 76, 024609 (2007). 
[28] A. S. Umar and V. E. Oberacker, Phys. Rev. C 73, 

054607 (2006). 



11 



[29] E. Chabanat, P. Bonche, P. Haensel, J. Meyer and R. 
Schaeffer, Nucl. Phys. A635, 231 (1998); A643, 441(E) 
(1998). 

[30] M. Kortelainen, J. McDonnell, W. Nazarewicz, 

P.-G. Reinhard, J. Sarich, N. Schunck, M. V. 

Stoitsov, and S. M. Wild, Ph ys. Rev. C (in press) 

flhttp: / /arxiv.org/abs/1 1 1 1 .4344 \ . 
[31] M. Kortelainen, T. Lesinski, J. More, W. Nazarewicz, J. 

Sarich, N. Schunck, M. V. Stoitsov, and S. Wild, Phys. 

Rev. C 82, 024313 (2010). 



[32] M. Kortelainen, J. Dobaczewski, K.Mizuyama, and J. 

Toivanen, Phys. Rev. C 77, 064307 (2008). 
[33] S. Landowne, C. H. Dasso, R. A. Broglia, and G. Pol- 

larolo, Phys. Rev. C 31, 1047 (1985). 
[34] H. Esbensen, S. H. Fricke, and S. Landowne, Phys. Rev. 

C 40, 2046 (1989). 
[35] K. T. R. Davies, K. R. S. Devi, S. E. Koonin, and M. R. 

Strayer, in Treatise on heavy ion Science, edited by D. 

A. Bromley, (Plenum, New York, 1985), Vol.3, page 3. 
[36] H. Esbensen, Phys. Rev. C 77, 054608 (2008). 



